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Introduction: 





Of all the species that inhabited the earth, humans shined with 
their uniqueness and their ability to combine and gather knowledge from their surroundings in 
order to lead a better life. Shelter or as we call 1t in the modern era “a building/a house”, was a 
necessity for the survival of our species. We needed a way to protect ourselves from the most 
destructive forces on earth, earthquakes, volcanos, tornados and many more. It was a hard task 
but throughout history, we made an impressive progress in building breath taking structures that 
withstand even the harshest climates. 


Structural design has always been a very interesting and creative segment in a large variety of 
engineering projects. Structures of course should be built in a way to resist applied forces stress 
and constraints and do not exceed certain deformations, in order to ensure safety. Moreover, 
structures should be economical. So certainly, the best designs are the ones that satisfy stress 
and displacement constraints and result in the least cost of construction. However, it's 
theoretical, and such designs only appear to exist on paper because we need to consider other 
factors. Consequently, our goal became the optimization of the weight, and the cost of 
structures while simultaneously maintaining safety when facing unpredictable factors such as 
natural disasters. 


Now if we dive deep in the nature of these forces, we find that they are dynamic with a 
magnitude that varies 1n time. Therefore, a dynamic analysis requires large amounts of data 
interpretation and information processing, resulting in a costly structure design. In a competitive 
world, optimum design has a great importance in an economic design of structures. 


* Ап optimizer is a formal plan or an algorithm, used to search for a “best” design that satisfies 
all conditions. 
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Choice justification: 


According to the national earthquake centre, there's an average of 20000 earthquakes yearly around the 
globe or approximately 55 earthquakes per day with 16 major earthquakes in the magnitude 7 range or 
greater. And since they are unpredictable, they can cause serious damage out of nowhere. Now this damage 
isn't caused by the quake itself but by the collapse of buildings with people inside them which is followed 
with deadly secondary effects and catastrophic aftershocks. Hence, earthquake proof buildings are a must, 
especially when no particular zone 1s safe from this natural occurrence. Seismic structural designs are based 
on a variety of interesting ideas that not only could save lives but also cut down on the costs of damaged 
buildings. 


The following picture 1s a map indicating earthquake zones with their intensity ranging from the least to 


most dangerous zones. 
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A general idea on how earthquakes affect buildings: 


Before we start optimising our designs, it is important that we first understand how earthquakes affect 
man-made structures. When an earthquake occurs, it sends shockwaves all throughout the ground, 1п all 
directions and within small intervals. Most buildings are generally made to resist vertical forces such as 
gravity and their weight, however they are poorly designed when it comes to side-to-side forces causing an 
almost immediate collapse to the structure when hit by strong seismic waves. Now you might be thinking 
that a skyscraper is more prone to collapsing, but it isn’t the case as tall buildings are more flexible and 
therefore can bounce back easier. 


When the ground shakes, buildings respond to the 
accelerations transmitted from the ground through the 
structure's foundation. The inertia of the building (it 
wants to stay at rest) can cause shearing of the structure 
which can concentrate stresses on the weak walls or 
joints in the structure resulting in failure or perhaps 
total collapse. The type of shaking and the frequency of 
shaking depends on the structure. 





How to build an earthquake proof building: 
Little story: 


When the shaking started at 5:46 a.m., Yasuhisa Itakura, an architect at a big Japanese construction company in 
Kobe, was sitting at his desk finishing a report he had toiled over all night. His office swayed, but the books 
stayed on their shelves and nothing fell off his desk. 


“T thought to myself, this earthquake 1s not that big,” Mr. Itakura said. 


It was, 1п fact, catastrophic. The Great Hanshin earthquake of January 17, 1995, killed more than 6,000 people in 
and around the industrial port city. 


Mr. Itakura had been cushioned from the violence of the earthquake because his three-story office building was 
sitting on an experimental foundation made from rubber — an early version of an engineering technique called 
base isolation. 


The technique that protected Mr. Itakura’s building is used in roughly 9,000 structures in Japan today, up from 
just two dozen at the time of the Kobe earthquake. Thousands of other buildings in the country have been fitted 
with shock-absorbing devices that can greatly reduce damage and prevent collapse. 


Chile, China, Italy, Mexico, Peru, Turkey and other countries vulnerable to earthquakes have adopted the 
technologies to varying degrees. But with notable exceptions, including Apple's new headquarters in Silicon 
Valley, the innovations have been used only sparingly in the United States. Seismic safety advocates describe this 
as a missed opportunity to save billions of dollars in reconstruction costs after the inevitable Big One strikes. 


To make a building more resistant to earthquakes, there are many methods: 


1. Create a Flexible Foundation 
One way to resist ground forces is to “lift” the building’s foundation above the earth. Base isolation involves 
constructing a building on top of flexible pads made of steel, rubber, and lead. When the base moves during the 
earthquake, the isolators vibrate while the structure itself remains steady. This effectively helps to absorb seismic 
waves and prevent them from traveling through a building, 





2. Counter Forces with Dampin 





You might be aware that cars have shock absorbers. However, you might not know that engineers also use them 
for making earthquake-resistant buildings. Similar to their use in cars, shock absorbers reduce the magnitude of 
shockwaves and help buildings slow down. This is accomplished in two ways: vibrational control devices and 
pendulum dampers. 


The effect of viscous damping and yielding, on the reduction of the seismic responses of steel buildings modelled 
as three-dimensional (3D) complex multidegree of freedom (MDOF) systems, is studied. The reduction produced 
by damping may be larger or smaller than that of yielding. This reduction can significantly vary from one 
structural idealization to another and is smaller for global than for local response parameters, which in turn 
depends on the particular local response parameter. 


The uncertainty in the estimation is significantly 
larger for local response parameter and decreases 
as damping increases. The results show the 
limitations of the commonly used static 
equivalent lateral force procedure where local and 
global response parameters are reduced in the 
same proportion. It is concluded that estimating 
the effect of damping and yielding on the seismic 
response of steel buildings by using simplified 
models may be a very crude approximation. 
Moreover, the effect of yielding should be 
explicitly calculated by using complex 3D MDOF 
models instead of estimating it 1n terms of 
equivalent viscous damping. The findings of this 
paper are for the particular models used in the 
study. Much more research is needed to reach 





more general conclusions. 


To withstand collapse, buildings need to redistribute the forces that travel through them during a seismic event. 
Shear walls, cross braces, diaphragms, and moment-resisting frames are central to reinforcing a building. 


Shear walls are a useful building technology that helps to transfer earthquake forces. Made of panels, these walls 
help a building keep its shape during movement. Shear walls are often supported by diagonal cross braces. These 
steel beams have the ability to support compression and tension, which helps to counteract the pressure and push 
forces back to the foundation. 


Diaphragms are a central part of a building's structure. Consisting of the floors of the building, the roof, and the 
decks placed over them, diaphragms help remove tension from the floor and push force to the vertical structures of 
the building. Moment-resisting frames provide more flexibility in a building's design. This structure is placed 
among the joints of the building and allows for the columns and beams to bend while the joints remain rigid. Thus, 
the building is able to resist the larger forces of an earthquake while allowing designers more freedom to arrange 


building elements. 
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In the following paragraphs, I'm going to explain the optimization method used for “reinforcing the building's 
structure". 


Definition: 


Steel moment framework uses a series of columns and beams, where the attachments are formed using a 
combination of welding and bolting. These connections are called "moment connections", with no releases to the 
joint when it is structurally analysed. The various connection methods used to achieve these connections have 
evolved throughout the decades, but as welding techniques and bolt strength have improved, the welded flange 
connection became the moment frame joint connection of choice. 


Optimisation Method: 


Eigenfrequency optimization describes a problem where one or more natural frequencies, of a structure or 
mechanism are optimized. This approach makes sense when the structures are subjected to dynamic loads whose 
response depends on one or more of its natural frequencies. Several problems involving dynamic loads in 
structural engineering have been solved using eigenfrequency optimization, such as maximizing the strength of 
RC column joints and vibration reduction in trusses. Additionally, this approach was used to optimize column 
dimensions using a surrogate membrane model of reinforced concrete moment frames, increasing the structure 
overstrength and ductility. 


The performance of a building during an earthquake seems to be greatly affected by the smallest natural 
frequencies of the structure. The reliance of the seismic response on the natural frequencies has been leveraged to 
develop efficient structural analysis procedures for calculating the seismic response of buildings, such as the 
modal pushover analysis. Moreover, several standards use the dominating nature of the smallest natural 
frequencies to recommend seismic response methods such as the equivalent lateral force (ELF) and the modal 
response spectrum analysis, both using elastic models of the buildings. 


The natural frequencies of a structure are properties that can be optimized using eigenfrequency optimization, 
suggesting that it is a suitable approach to design buildings with improved seismic performance. Furthermore, the 
optimization of the natural frequencies of a structure has low computational costs since the objective function uses 
the elastic model of the building to make the calculation, which is significantly less complex than the nonlinear 
analyses required by existing methods for optimizing the design of SMRF. 


Developing an efficient formulation of eigenfrequency optimization for SMRF has three challenges. The first 
one is that the elements sections in a SMRF are represented by discrete variables, given by the typical steel section 
profiles used in construction practice. A second challenge stems from the fact that the problem optimizes natural 
frequencies; therefore, it is necessary to introduce a constraint to the problem that limits the cost of the building, 
ensuring a competitive design. A third challenge is ensuring that the resulting design satisfies design code 
requirements for resistance and functionality. 
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Handling this problem requires expressing the element sections of the SMRF as a vector of integer indexes, 
each one corresponding to a section profile, which has associated the corresponding section properties (area, 





Challengel: 


moment of inertia, etc.) 
Challenge 2: 


This is addressed by imposing a limit WO to the total weight of the building W(p), which is a function of the 
profiles p of the structural elements. 


Challenge 3: 


To overcome this challenge, each element of the building (1.e., columns and beams) must be limited between 
a minimum profile (Pi; min) and a maximum profile (Pi; max). In addition to fulfilling design code 
requirements, this scheme also contributes to limit the search space, avoiding the possibility of generating 
unviable solutions that would be discarded by the design code constraints. 


With those considerations into account, the equation below shows the formulation of the eigenfrequency 
optimization problem for SMRF, based on the stiffness (K) and mass matrixes (M) of a given structure. Here 
Wn, is defined as the n-th natural frequency of the building and P is an integer vector with the section profiles, 
which has associated the corresponding section properties: 


max ШАРЕ, 
subj ect to Pimin = р = Pimax 
W (ip) < Wo 


[К — cM | Q,-0 


This formulation seeks to maximize the set of N first natural frequencies of the structure, subjected to three 
constraints. In the first constraint, Pi; min and Pi; max are the minimum and maximum index values of each 
element represented in P. The second constraint accounts for the volume limit previously discussed. The third 
constraint comes from the structural dynamic theory, and it means that Wn 15 a natural frequency of the 
structure. 


Evaluation of Set of member 
gravity and dimensions is 
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Feasibility, understood as the degree to which the optimization procedure is easy or convenient to implement, 
is one critical aspect of an optimization method to succeed in engineering practice. As demonstrated in equation 
above, the proposed method uses the elastic properties of the building, so it 1s feasible for implementation in 
engineering practice. 


This approach has several advantages. First, the design that results from the optimization satisfies not only the 
strength requirements but also the drift limit, as the method is maximizing the fundamental frequency (1.e., 
minimizing the fundamental period) and the drift limit is checked in the elastic range. Consequently, optimized 
building will have a smaller one than the traditionally designed building. The second advantage is that unlike 
existing methods, the proposed method does not need to perform any nonlinear analyses, significantly improving 
the computational performance. A third advantage of the proposed method is that the proposed formulation 
works for any type of buildings, regardless of their irregularity and whether they are two-dimensional and three- 
dimensional problems. In the latter case, the stiffness matrix must be formulated using the frame-type element 
which has 6 degrees of freedom on each element end. Similarly, the mass matrix must account for the additional 
degrees of freedom. In all cases, the method requires that the strength and ductility requirements must be met by 
the proper selection of the Pi; min and P1; max the for the members profiles. 


Figure 1 :a four-story building with three identical beam spans of 6 m, and all the floors have a story height of 3 m. 


This building was designed per the Load and Resistance Factor Design in the AISC 360-16 following the 
customary engineering practice. Our goal is to optimize the design of this structure following the method stated 
above. 


Two section profile databases were used for this purpose. Tables 1 and 2 show the profiles used for beams and 
columns, which range from the IPE100 to the IPE600 and from HEA100 to HEA600, respectively. 
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Figure 2': IPE100 to IPE600 


2,90 429 11.2 420 62.2 3.02 


5.360 anf 12.5 604 50.5 3.35 





Profile ID Inertia 
(strong 


axis) cm? 


Inertia 
(strong 


Inertia Profile 


(strong 


Profile 


axis) axis) cm? 


cm? 

































































HEA100 1 21.2 | 349.2 HEA280 10 97.3 136699  HEA600 19 226.5 141199.9 
HEA120 25.3 606.0  НЕАЗОО (11 112.5 182601 
Figure 2" : HEA100 to HEA600 
 HEA140 - 31.4 1033.1 НЕАЗ20 12 1244 . 229302 
HEA160 38.8 | 1672.8 |  HEA340 13 133.5 27689.8 
HEA180 45.3 2509.9 HEA360 |14 142.8 33090.0 
 HEA200 53.8 | 3692.0 HEA400 15 159.0 45070.0 
` HEA220 (643 5410.2 HEA450 16 178.0 637200 
HEA240 76.8 | 7763.1 | HEASOO 17 197.5 86969.9 
 HEA260 > 86.8 104499 HEAS50 18 211.8 1119001 


The objective function and constraints were coded in MATLAB® and the multiobjective evolutionary 
algorithm, Borg-MOEA was used to solve the problem. 


The proposed method is applied to a four-story building. This building ıs first designed according to the 
customary engineering practice for a base shear of 161.5 kN. The design results show that using profiles ID = 11 
(HEA300) and ID = 10 (IPE100) for columns and beams satisfy all the drift, strength, and deflexion 
requirements. This design has a mass of 7.28 ton, which is used as WO the limit for the problem formulation. 
Since this 1s a low-rise building, the fundamental frequency was considered as the only objective function. 


A key step to the optimization of this structure is choosing suitable values for Pi; min and Pi; max and defining 
the first constraint. 


Pi; min: a good strategy for beams is selecting the smallest profile that satisfies the serviceability (deflexion) 
requirements of the design code. One way to conduct this task 1s analysing the floor system for the gravity loads 
and finding the profile for which the deflexions are closest to the code limit for the building. For this structure, 
beam ID = 9 (IPE270) satisfies the strength and deflexion requirements; hence, they are used as the Pi; min. 


We should note that Pi; min should be selected simultaneously with Pi; max, and for this case column ID = 10 
(HEA280) and beam ID = 11 (IPE330) are selected for the columns Pi; min and beams Pi; max. 


The table below shows the results of the building design by the traditional method and the optimizing method: 


Traditional desien Optimized desien 
Story eam 
Column ID Beam ID Column ID ID 
1 11 10 12 10 
2 11 10 12 10 
3 11 10 10 10 
4 11 10 10 10 


Conclusions: 


This table shows that the proposed method did not change the building's beam configuration, leaving all floors 
with the initial ID = 10 (IPE300) profile. This result may stem from the additional constraint imposed to the 
objective function that beams within a same floor must use the same section profile. This constraint was 
included to accommodate construction practices, where it is common for all beams within a same floor to use 
the same section. Regarding columns, the method assigns ID = 12 (HEA320) sections to columns of the first two 
stories and ID = 10 (HEA280) to columns in the top half of the building. The material consumption for this 
building is 7.22 ton, a 0.8% reduction compared to the traditional alternative. In terms of ratio of the column 
inertia to beam inertia, column ID = 12 to beam ID = 10 in the first two floors of the optimized design is 2.74, a 
25% increase compared to the 2.18 in the traditional design. On the other hand, for floors 3 and 4, this ratio 
decreases to 1.63. In all cases, these ratios satisfy the code constraints. 


The results for the traditionally designed building and the building designed per the proposed method аге 
shown in the picture below. The X-axis shows the roof drift ratio (КОК) of the buildings as percent, and the Y- 
axis shows the ratio (Vs/Vdesign) between the base shear Vs 1n the pushover and the design base shear 


Vdesign of the building. 
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Compared to the traditionally designed building, 
the seismic performance at the life safety and 
collapse levels for the building designed 
following the procedure presented is improved. 
These benefits come without any additional 
costs in materials and were achieved after less 
than 30 minutes of computational time. 


The results from this analysis show that the 
buildings have almost identical performance up 
to 1.5% roof drift, suggesting that the seismic 
performance for small earthquakes will be 
similar. 


After this point, the building designed per this 
method shows better performance, exhibiting a 
higher overstrength. 
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Compared to existing alternatives, this method based on eigenfrequency optimization has the 
advantages of being computationally efficient and feasible to implement. These benefits stem from the 
fact that the method uses the information of elastic models of the buildings. 


This method produces designs of steel moment frames that have superior seismic performance than 
buildings designed following the customary engineering practice. The results from the seismic 
performance for a four-story building show that the building designed using the proposed method is 
12% less susceptible to collapse than its traditionally designed counterpart. Eigenfrequency 
optimization 1s a viable framework to optimize the design of steel moment frames, which lends itself 


to be implemented and solved using evolutionary approaches like genetic algorithms. The method can 
be used for buildings with more complex dynamic responses as it allows incorporation of more than 


one natural frequency as objective function. 


Perspectives of expanding the subject: 


The seismic activity in Morocco is not as important as in other world areas such as Japan, California, Greece, 
Turkey, etc. Nevertheless, background seismicity in Morocco is not negligible and has produced a few 
historical and instrumental earthquakes with local magnitudes above 6. The Agadir earthquake 1960, for 
instance, destroyed this city with several villages and caused 12,000 deaths. Later, other remarkable 
earthquakes occurred: in Al Hoceima in 2004 in the North of Morocco, caused 629 deaths. However, seismic 
risk is relatively high due to vulnerable buildings that are either old or don't respect seismic standards. One of 
the earliest projects regarding earthquake proofing buildings in morocco date back to 1960 after the deadly 
earthquake in Agadir and it was called “building a new resistant city”; A High Commission for Reconstruction 
was established and special earthquake resistant methods were officially approved. No building permission is 
given if even the smallest construction does not conform to the “Norms Agadir 1960.” These norms define all 
construction requirements and stress the importance of avoiding unsafe building sites such as steep slopes, 
water- impregnated ground, sandy soil or rubble filled soil. Rocky ground is considered best. The 
specifications call for horizontal foundations but if the foundation is built on two levels, the joints should be 
reinforced. There must be steel wire fabric reinforcement at ground level. Joints must also be provided be-t 
ween the two sections of a building for they will react differently to seismic vibrations. Even the type and 
quality of building materials such as sand, bricks, steel, cement and water used in mixing concrete are defined. 


Since the concept isn't new but certainly needs improvements, the Eigenfrequency optimization is the perfect 
way to both optimize the weight and the cost of future seismic designs. We can start by proofing small 
portions of zones of the highest risks and this means Seismic retrofitting old houses and optimizing their 
weights while taking in consideration forces of stress, constraints and side to side forces. 


Future research on this topic can be pursued in several topics. One interesting area is evaluating the benefits of the 
proposed method for other building configurations, particularly for high-rise buildings where several frequencies are 
needed to capture the dynamic response and which would be required as objective functions. Considering the 
potential to be a computationally efficient method, ії is a worthy topic to investigate the method application to 
optimize three-dimensional buildings with different types of irregularities. 


Conclusion: 


This study employs an organized methodology for finding the optimum shape of the structures under earthquake 
loading. The numerical examples show that the optimum shapes are obtained in a few iterations. Furthermore, it could 
be applied easily with small modifications to three-dimensional structures. To finish this research, | would like to add 
that this is not the only method available to optimize the shape of the structure under earthquake or generally 
speaking dynamic loads. As mentioned before, there are many ways to build a resistant structure and therefore there 
are many algorithms to do so. | would like to mention the method of the finite element analysis where earthquake 
loading is transformed into the equivalent static loads (ESLs) and the stresses are computed from the differentiations 
of the displacement field. 


As technology goes further, we may be able to create different and more efficient ways to optimize our structure and 
probably even predict the time of the occurrence of an earthquake leading to the creation of better shields that 
protect future generations from our furious earth. 
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